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Abstract 

We derive a solvable resonant-level type model, to describe an impurity spin 
coupled to zero-energy bound states localized at the edge of a one dimensional 
d-wave superconductor. This results in a two-channel Kondo effect with a 
quite unusual low-temperature thermodynamics. For instance, the local im- 
purity susceptibility yields a finite maximum at zero temperature (but no 
ln-divergence) due to the splitting of the impurity in two Majorana fermions. 
Moreover, we make comparisons with the Kondo effect occurring in a two 
dimensional d-wave superconductor. 

PACS. 71.10.Pm - Fermions in reduced dimensions (anyons, composite 
fermions, Luttinger liquid, etc.). 

PACS. 74.20.Mn - Nonconventional mechanisms (spin fluctuations, polarons 
and bipolarons, resonating valence bond model, anyon mechanism, marginal 
Fermi liquid, Luttinger liquid, etc.). 

PACS. 72.15.Qm - Scattering mechanisms and Kondo effect. 
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Introduction and motivation. — The problem of magnetic impurities in a d-wave super- 
conductor has been extensively discussed in the literature. This problem is of relevance for 
high-T c and heavy-fermion superconductors. Since the pioneering paper by Abrikosov and 
Gor'kov m it is known that magnetic impurities act as pair-breaking centers 0|| and hence 
reduce the amplitude of the condensate. Recently, it has been shown that quantum me- 
chanical fluctuations of magnetic impurities give rise to other important effects, especially 
when coupled to the fermionic quasiparticles of the superconducting state. In a conventional 
metal, these correlations lead to the so-called Kondo effect [|J, i.e. the complete screening 
of the impurity at low temperature. In d-wave superconductors, the density of normal 
quasiparticles vanishes at the Fermi level. Then, screening is absent in perturbation theory 
and a critical exchange coupling between the quasiparticles and the magnetic impurity is 
crucial for Kondo screening to take place 0. Using an effective model of one-dimensional 
chiral fermions coupled to a single impurity and the so-called large N-expansion (N, being 
the rank of the symmetry group of the impurity spin), Cassanello and Fradkin have shown 
that there is a quantum phase transition || from an unscreened impurity state to an over- 
screened Kondo regime [[?[] . This occurs at a critical value J c of the coupling constant which 
is typically a fraction of A, the size of the gap away from the nodes at zero temperature. 
In the overscreened phase, which arises due to the consequent coupling between normal 
quasiparticles close to one particular node on the Fermi surface and the impurity, the low- 
temperature thermodynamics is given by Ximp ~ TlnT and Cj mp ~ T 2 lnT. In this regime 
the Kondo scale Tf. is almost always smaller than A, such that pair-breaking effects should 
not be prominent ||. 

On the other hand, the spin liquid phase of the two-leg ladder system gives us a sim- 
ple model of a Mott-insulating phase which actually exhibits pairing, with an approximate 
d-wave symmetry. Moreover, upon doping, the two-leg ladder yields quasi-long-range su- 
perconducting (d-wave) pairing correlations [§,|K|]. Experimentally, superconductivity has 
up to now only been observed in a two- leg ladder material under high pressure In this 
Letter, we investigate for the first time the (two-channel) Kondo model occurring in a two- 
leg ladder, for a weak repulsive electron-electron interaction U. First, based on bosonization 
methods and previous results, we derive an effective Hamiltonian (called below, two-band 
Luther-Emery Hamiltonian) describing well a two-leg ladder d-wave superconductor at zero 
temperature. Second, we show that the Emery-Kivelson solution of the usual two-channel 
Kondo model [|l^] and the effective two-band Luther-Emery Hamiltonian, can be combined 



leading order in U. The model becomes solvable and the low-temperature thermodynamics 
can be computed explicitly. Additionally, an exact solvable model for an impurity coupled 
to a one- dimensional s-wave superconductor has been proposed in ref. In comparison, 



the problem of a magnetic impurity in Luttinger liquids has been more studied |14|jn 

Model. — To study low-energy physics, we use a continuum limit formalism (a is a required 
short- distance cutoff). The noninteracting two- leg ladder of spinful fermions is described by 
the Hamiltonian density 

Ho = -tJ2tf ja (x + a)V ja (x) + h.c. -t ± J2^\M^-ia(x) + h.c, (1) 

where t ant t± are the hopping amplitudes along and between the chains and tyj a {x) creates 
a fermion with spin a =1,1 in the chain j at the rung x. In what follows, we are going 
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to consider the limit U <C t and t± ~ t. Therefore, we can diagonalize first H in the so- 
called bonding-antibonding band basis, defining: ^ oajWa = l/v2( 1 I r i a ,± 1 i r _ la ). The bonding 
\l/ oa and the antibonding ^ na operators are expanded as (j = o,ir): \F, a ~ ^Rj a e lkFjX + 
Lja£~ lkFjX ■ Upon linearizing the spectrum around the four Fermi points the kinetic energy 
takes the form, 

Ko = -v £ l^ Rja id x q Rja -V{ ja id x y Lja ]. (2) 

j=o,7r;a 

Moreover, not so far from half-filling the Fermi velocity in each band is the same, hereafter 
denoted as v Now, we must include the effects of the electron-electron interactions. 

As a generalization of the Hubbard interaction, the two-leg ladder is described by a set of 
marginal Fermion interactions (consult Eq.(2.5) of ref. []T0|] ) , that can be readily decomposed 
in terms of the continuum Dirac fields. Remarkably, the renormalization group equations 
reveal that before the flows of the different parameters scale towards strong coupling, they 
are attracted to a special direction, characterized by a unique coupling constant g ~ U . 
Low-energy physics depends only on this coupling. To describe it precisely we need to 
"bosonize" the electron fields as follows Jl7| 



*p ja = « ia e^- +e -\ (3) 

with P = R/L = ± and Kj a the corresponding Klein factors. The displacement field <f)j a 
and phase field Qj a satisfy the usual commutation rules. It is convenient to separate charge 
and spin modes, defining: <f>j C = (0^ + 4>j\) / V% and 4>j S = — (f)ji) / \/2. The Hamiltonian 
density reads (see Eqs.(3.2) and (3.4) of Ref. |T§): H = H + Hi, with 

n o = \ E [^{dAaf + vK{d x Q a )% (4) 

a=jc,js 

Hi = —g cos v^2(6 oc — Q nc ) cos \/2(f) os cos y/2<f) ns , 

K = 1 — g/nv < 1, being the Luttinger parameter of the effective theory. Hi is produced 
chiefly by interband charge and spin Cooper backward scatterings. 

When the single coupling constant g is flowing off to strong couplings, the two spin fields 
cf) os and <f) ns will be pinned in the minima of cosine potentials, i.e. <ftj S ~ 0. The spin in each 
band vanishes |18| . Furthermore, the coherence between the bands arises due to the pinning 
of Q c f = (O oc — Q^c) j \pl ~ that guarantees the binding between 'holons' of the two bands 
the superfluid system is ruled by the unique phase C = (0 OC + Q nc )/y/2 = \/2<S)j C . As 
a consequence, the pair field operator Aj = ^Rj^Lji ~ ±expi6 c yields only a very slowly 
decreasing power-law decay: < A^(a;)Aj(0) > ~ x~ l ^ 2K . This produces prominent d-wave 
superconductivity ||20|| . Since K < 1, the (marginal) operator in Hi is relevant and then 



generates a superconducting crossover temperature A = t exp(— nv/g) (or the spin gap at 
the Fermi level) [1B1. In this Letter, we decouple Hi = Hb s + Hhoion, with: 



H bs oc -#(cos V8(p os (x) + cos V8(f) ns (x)), (5) 
H ho i on oc g cos y/8Q c f (x) . 
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This decoupling does not affect the ground state properties and the gap equation. Moreover, 



the effective spin terms in Eq.(y) produce explicitly a two-band Luther- Emery model pi 
The model we propose here to describe is an impurity spin situated at the edge x = of such 
d-wave superconductor. The bulk Hamiltonian obeys: H = / L dx 7i. In the thermodynamic 
limit we have a semi-infinite system. 

First, it is convenient to rewrite the spin part H os + Hb s on the entire axis — oo < x < oo, 
redefining ^Rj a {x) = ^fj a (x) and ^Lj a (x) = ^j a (—x), with i&j a (x) ~ exp i<&j a (x). In terms 
of bosonic operators, one can simply equate (see Eq.(|3|)): 

(<f> js (x), Q js (x)) = ^ js (x) =F ^js(-x)). (6) 



Rescaling the variables in a standard manner <pj S — > V K<pj S and Qj S — > Qj S / V K, the (bulk) 
spin Hamiltonian becomes: 

y /*+00 

H os = - dx (d x $ os ) 2 + (d x $ ns ) 2 (7) 
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oo 

+oo 



U, — y / rl-r p iV2K<5> 03 (x) -iV2K<5> os (-x) -iTrKsgn(x) 

bs ~ (27ra) 2 i-oo +l j ' 



For a complete introduction on bosonization with boundaries, consult for instance ref. |22 
Second, in the band basis the Kondo interaction reads 

H k = J - £ lim_ {*y*)a^(x)} • 5 = J J2 lim x ^ 0+ {^ Rja (x)a^ RjP (x)} ■ S. (8) 

Backward scattering processes are already implemented by the open boundary: [^ija{0) + 
^Rja(0)] = jl3fl . Using preceding transformations, the electron fields of Eq.(||) read: 



Hf Pja ( x ) = 9 3a (Px) = « iaC «( pa ^(»i-(*)-*i-(-*)Hi^(»i«W+*i«(--))). 
The vanishing of the wave function at the origin gives also: 



(9) 



and 22 



lim x _ 0+ {$ jQ (x) - ^ja(-x)} = ti/Vk (10) 
\im x ^ 0+ {$ ja (x) + $ ja (-x)} = 2$ iQ (0) - tv/t/K, 



^i?ja(0) = lim x ^ +^ ia (x) ~ K ia e . (11) 
Therefore, the two-channel Kondo interaction takes the form [ 23[| : 

= ^=^6^(0) + — cos $ sf (0) exp(z$ s (0))5 + + £>(— ) + h.c. (12) 

2nvK ira irv 



We have canonically transformed $ s = ($ os + ^^ s )/y/2 and $ s/ = ($ os - $ ws )/y/2. The 
charge degrees of freedom are not coupled to the impurity, so we subsequently omit them 
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in the following. Moreover, it is convenient to use the unitary rotation: U = exp(i$ s (0)S' 2 ). 
We can first absorb the z-component at the particular point J z = -KVyK, namely the 
generalized Emery-Kivelson point. Then, the planar Kondo term is transformed as in the 
usual two-channel Kondo model |i"2|j24fl : 



J 



H k = UH k U ] = — cos $ s/ (0) • S x + h.c. 
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Leading order in g, the backscattering Hamiltonian becomes: 

H bs = H bs oc +g f + °°dx e <*.W e -<*.(-*)[ e »./(*) e -<*./H + e -»./(*) e »./(- 

J — oo 



xh 



(13) 



(14) 



Finally, we must retransform the bosonic variables as spinless fermions. Using the Fermi 
operators ~ expi$ s (x) and ^ s f ~ expz$ s /(x) fi~2| , |24| 



, one arrives at 



r+oo 

H bs = ~g dx Vl(x)y s (-x)yl f (-x)y sf (x), 

J -co 
r+oo 

— —g I dx ic ^l(x)^ s (—x)sgn(x) + (s — > sf). 
Using ref. pjf] , one definitely obtains H = H os + H bs + H k , with 



(15) 



H os + Hi 



6.s 



E 

j=s,sf 



dx [v^](x)(-id x )^j(x) - i&$} j (x)fy j (-x)sgn(x)] 



(16) 



Hu 



J 



[* sf (0)+*l f (0)]a + h.c. 



(17) 



S x = a/y/2, S y = i>/\/2 and S z = tab. The a and b are Majorana (real) fermions. Thus, 
we arrive at a resonant-level type model that is solvable because it combines the Emery- 
Kivelson solution of the two-channel Kondo model with the two band Luther-Emery model. 



Now we introduce another important energy scale, the bare Kondo resonance |T2| , |24| , |26 
T = J 2 /irva. It is instructive to diagonalize first the bulk Hamiltonian, similarly as in ref. 
|i~3| . Bulk solutions for quasiparticles are plane waves exp ±ipx, where e p = ±[(vp) 2 + A 2 ] 1 / 2 . 
But, starting with a semi-infinite system should produce novel Kondo features. Precisely, 
we have new solutions at the edge x = ('bound S=l/2 spinons') with e p = i.e. p = 
±iA/v. Conversely to usual superconductors these zero-energy bound states exist without 
any proximity effect with a normal metal [p7]. The impurity Green's functions take the 
forms: 



1 p 
A(r) = -(Ta(r)a(0)) = -E 



T 



1 e 
■(Tb(r)b(O)) = -£ 



iWnT 



lW r 



A(iw n ) 



(18) 



w n are fermionic Matsubara frequencies and A(iw n ) = T(y/w n 2 + A 2 + A)/iw n 
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Physics and conclusions. — A new interesting physics arises when V <C 2 A. The normal 
quasiparticles definitely disappear from the spectrum at low energy pQ |. A zero-energy 



spinon, from the s /-sector, should form a stable bound state with half of the impurity spin: 
the associated (rescaled) Kondo energy scale is w co h = (2TA — T 2 ) 1 / 2 , solution of iw n = 
A(iw n ). The situation becomes completely different from that in a two-dimensional d-wave 
superconductor where quasiparticles near the nodes mainly screen the impurity. In that case, 
both Ximp and Ci mp go to zero at T = meaning that the impurity is completely screened. 
Furthermore, the reason behind the existence of J c is that the effective interaction between 
the impurity and the normal excitations is momentum dependent and that it vanishes at 
small momenta ||. In our case, an overscreened Kondo effect can already take place when 
J — > 0. When TcT, integrating over the coupling constant the free energy takes the form: 

F = F Q + F(a) = F - | In cosh (19) 

For J = 0, we have S = —dF /dT = In 2 because an isolated impurity has a two fold 
degenerate ground state. Then, expanding the second term of Eq.(|T9"[) near T — > 0, the 
result is S(a) = — In 2/2. We find the same residual entropy S(0) = S Q + S(a) = In 2/2 as in 



the usual two-channel Kondo model (in a metal) |L2[p4 |: only half of the impurity is screened 
at zero temperature. Furthermore, the impurity contribution to the specific heat comes from 
thermal activation, i.e. C imp oc (1/T 2 ) • exp — w co h/T Ul3 |. Another important point is that 



performing a deviation from the generalized Emery-Kivelson limit does not change much 
this conclusion: the Green's function for a bound spinon is of the form of B(t), and then 
no ln-divergence may arise in the second order correction (in A = J z — tcv\^K) to the free 
energy. We can also compute the local impurity susceptibility (i.e. when the magnetic field 
acts only on the impurity spin) ||28|| . When T <C T, it is convenient to use the following 
asymptotic form of the propagator (0 < r C 

A(t) = - e ~ w °° hT . (20) 

This leads to 



' i 



Xi- 



imp {T) = fair (S z (t)S z (0)} cx — { e —^ _ e -WT } . (21) 

J To w coh 



The exponential contained in a correlations cuts the r-integral off (r G ~ 1/T). The most 
interesting observation is the absence of a ln-divergence in Xi-imp a t zero temperature, despite 
the splitting of the impurity into two Majorana fermions a and b |28|| . It is noteworthy that 
the splitting phenomenon produces only a finite maximum in the local impurity susceptibility 
when T — > 0. When T ~ 2A, the bound state formed between the Majorana fermion a and 
an edge spinon disappears i.e. w co h = and for larger couplings only normal quasiparticles 
will be involved in the Kondo coupling with the impurity. At low temperatures, this results 
in C imp {T) and xi-imp(T) 1/T (obtained by expanding Eq . (|2l|) near w co h = with 
T o — * 0). We have a crossover region where the impurity is completely unscreened; the 
density of states at the Fermi level tends to zero and quasiparticles cannot pretend to form 
a stable bound state with the impurity. Finally, in the limit A — >• (which corresponds to 
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K=l), the density of states at the zero-energy level becomes finite and we recover previous 
results: A(iw n ) = —iTsgn(w n ) and then Xi-imp(T) oc InT, C imp {T) oc TlnT [Q. The 
impurity is now overscreened by two channels of normal electrons. 

To conclude, a solvable model describing a magnetic impurity coupled to bound spinons 
located at the edge of a one dimensional d-wave superconductor, is proposed in this Letter. 
We find a new type of two-channel Kondo effect resulting in a quite unusual low-temperature 
thermodynamics. 

I thank the theoretical group of ETH for stimulating discussions, and especially U. 
Ledermann and M. Rice for collaboration on a related work. This paper is devoted to my 
father. 
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